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P ! ABSTRACT 

Context. It remains difficult to interpret the near-IR emission of young stellar populations. One main reason is our incomplete 
1 1 , understanding of the spectra of luminous red stars. 

■ Aims. This work provides a grid of theoretical spectra of red giant and supergiant stars, that extends through optical and near- 
' | IR wavelengths. For the first time, models are also provided with modified surface abundances of C, N and O, as a step towards 

■ accounting for the changes that occur due to convective dredge-up in red supergiants or may occur at earlier evolutionary stages 
^v^j in the case of rotation. The aims are (i) to assess how well current models reproduce observed spectra, in particular in the 

I , near-IR, (ii) to quantify the effects of the abundance changes on the spectra, and (iii) to determine how these changes affect 
^v^j ■ estimates of fundamental stellar parameters. 

Methods. Spectra are computed with the model atmosphere code PHOENIX and compared with a homogeneous set of observations. 
Although the empirical spectra have a resolution of only A/AA ~ 1000, we emphasize that models must be calculated at high 
spectral resolution in order to reproduce the shapes of line blends and molecular bands. 
| Results. Giant star spectra of class III can be fitted extremely well at solar metallicity down to ~3400 K, where difficulties appear 
• • . in the modelling of near-IR H2O and TiO absorption bands. Luminous giants of class II can be fitted well too, with modified 
. ' sur ^ ace abundances preferred in a minority of cases, possibly indicating mixing in excess of standard first dredge-up. Supergiant 
k> , stars show a larger variety of near-IR spectra, and good fits are currently obtained for about one third of the observations only. 
» I ■ Modified surface abundances help reproducing strong CN bands, but do not suffice to resolve the difficulties. The effect of the 
C\3 abundance changes on the estimated T c ff depends on the wavelength range of observation and can amount several 100 K. 

Conclusions. While theoretical spectra for giant stars are becoming very satisfactory, red supergiants require further work. The 
model grid must be extended, in particular to larger micro-turbulent velocities. Some observed spectra may call for models with 
even lower gravities than explored here (and therefore probably stellar winds), and/or with more extreme abundances than 
predicted by standard non-rotating evolution models. Non-static atmospheres models should also be envisaged. 

Key words. Stars: fundamental parameters - Stars: red supergiants - Stars: red giants - Stars: atmospheres - Stars: spectra 



^1- 

o 



1. Introduction 

Red supergiants and red giants are the most luminous 
stars in, respectively, star forming or old passive galax- 
ies. Being cool, they are the dominant sources of near-IR 
light. In highly reddened starburst galaxies, the near-IR 



Send offprint requests to: A. Langon e-mail: 
lancon@astro . u-strasbg . f r 

* Selected theoretical spectra (see text) can be 
retrieved in FITS format at CDS via anonymous 
ftp to cdsarc.u-strasbg.fr (130.79.128.5), or via 
http://cdsweb. u-strasbg. fr/cgi-bin/qcatTJ/A+A/] 



light from red supergiants is sometimes the only direct 
information available on the stellar populations. Models 
for their spectra are thus important, even though they are 
also particularly difficult to construct (rich molecular line 
spectrum, extended atmospheres). If they are successful 
in reproducing empirical spectra, it will be legitimate to 
use them instead of observed spectral libraries in future 
analyses of galaxies. 

Recently, Levesque et al. (2005) have shown that up- 
to-date models compare well with optical observations 
of red supergiants, and shown that this success helps 
in explaining the location of observed red supergiants 
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in the HR-diagram. At near-IR wavelengths (1-2.5/im), 
the most prominent molecular features are those of CO 
and CN. Their sensitivity to surface gravity and effec- 
tive temperature (T eS ) has been the basis of the 8-colour 
classification system of White & Wing (1978), although 
they are also sensitive to other parameters (Tsuji 1976, 
McWilliam & Lambert 1984, McGregor 1987, Bessell et al. 
1991, Origlia et al. 1997). Strong CN bands are predicted 
for low gravity stars with temperatures around 4500 K, 
and are indeed observed in local red supergiants (Langon 
& Wood 2000) and in extragalactic objects such as the 
bright star clusters of M82 (Langon et al. 2003). However, 
it had not been verified until now whether models are ca- 
pable of reproducing the strengths of various CN and CO 
bands throughout the near-IR range simultaneously. Nor 
whether they can match optical and near-IR properties 
together. 

An important aspect not accounted for in recent collec- 
tions of model spectra for red supergiants is internal mix- 
ing. Standard stellar evolution predicts non-solar surface 
abundance ratios due to convective dredge-up in the red 
supergiant phase (Iben 1966, Maeder 1981). Early obser- 
vations had pointed out the inadequacy of solar abundance 
ratios in individual cases (e.g. a Ori, Beer et al. 1972). 
More recently, both theory and observations showed that 
main sequence rotation or other processes are capable 
of mixing CNO-cycle products into the atmosphere even 
before the red supergiant phase is reached (Maeder & 
Meynet 2001, Trundle & Lennon 2005). In red supergiants, 
He and 14 N surface abundances are typically enhanced 
while 16 and 12 C abundances are reduced. Modified 
abundances of C,N and O alter the relative strengths of 
the predominant molecules. 

In this paper, we present recent PHOENIX models specif- 
ically computed to address some of the above points. The 
emphasis is on the effects of non-solar abundance ratios 
of C,N and O; a more complete study of other param- 
eters (in particular micro-turbulent velocities) has been 
started and will be reported in a forthcoming paper. The 
model assumptions are described in Sect. [5] and the pre- 
dicted colours and molecular features in Sect. [3] In Sect. [4] 
and[5l the models are compared with spectroscopic data 
covering wavelengths from 1 to 2.4 /im or, for a subsam- 
ple, from 0.51 to 2.4 /im. Giants of class III, luminous gi- 
ants of class II and supergiants of class I are discussed 
successively. The discussion in Sect. [HI focuses on funda- 
mental parameter determinations from spectra, including 
the effects of mass and surface abundances on T eff . A brief 
conclusion summarizes the results. 

2. PHOENIX models with solar and modified 
abundances 

2.1. Summary of model ingredients 

The model atmospheres and synthetic spectra were com- 
puted with PHOENIX version 13.11.00B. The model setup 
is identical to that of Kucinskas et al. (2005). We recall 
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Fig. 1. Typical differences between PHOENIX spectra 
obtained with an initial wavelength sampling step 
of 2 A (grey) and 0.1 A (black). Both spectra have 
been smoothed by convolution with a gaussian with 
FWHM=15A. The models shown have T cff =4000K, 
log(g)=l, M=1M Q , but differences are important for any 
of the calculated models. Only the high resolution calcu- 
lations match the data. 



only the most relevant details here. The models are com- 
puted in spherical symmetry. A mixing length to pressure 
scale height ratio of 2.0 is used for all models. Dust is 
allowed to form in the atmospheres but is assumed to im- 
mediately rain out of the photospheric layers; therefore, 
no dust opacities are used in the models shown here. This 
is an important assumption for cool models with large ex- 
tensions. In addition, all models presented here have large 
enough gravities to not produce a radiatively driven wind 
and, therefore, winds are not included. 

The model spectra were computed specifically for com- 
parison with data that has a spectral resolving power of 
order 1000, i.e. AA ~ 10 A at 1 /im. Nevertheless, we em- 
phasize that the model spectra must be computed at high 
spectral resolution before smoothing, in order to sample 
individual absorption line profiles properly and to obtain 
low resolution spectra that resemble observations (Fig.Q]). 
We used a wavelength sampling step of 0.1 A throughout. 
Using only half these points produces negligible changes at 
A > 8000 A (0.1 % rms), and small variations in the shapes 
of the strongest optical bands at A < 8000 A (2% rms). 
The small sampling step used here is an important change 
with respect to previous collections of PHOENIX spectra, 
which were computed with an initial wavelength sampling 
step of 2 A (e.g. Kucinskas et al., 2005, and models in- 
cluded in the library of Martins et al., 2005). 
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The models discussed cover effective temperatures 
ranging from 2900 to 5900 K and gravities in the set 
log(<7)={— 1,-0.5,0,1,2} (cm.s -2 ). The micro-turbulent 
velocity is set to a constant value of v llllc — 2kms~ 1 , ex- 
cept in a few exploratory models. Values of 2 to 3kms _1 
are typical for red giant stars (Smith & Lambert 1985). 
A more extensive grid of models covering higher values of 
this parameter is in the process of being calculated. Two 
stellar masses are considered: 1 M© and 15 Mq. Models 
at 9 Mq were also computed, but the differences with the 
15 M Q ones are negligible. For M=1M Q , many of the cal- 
culations at log(<7)=— 1 did not converge (radiation pres- 
sure incompatible with the assumption of no winds), and 
this mass-gravity combination is therefore excluded from 
the discussion. We also restrict the discussion to optical 
and near-IR wavelengths, with a focus on wavelengths be- 
tween 0.81 and 2.4 /im. 

2.2. Abundances 

The reference set of models assumes solar abundances, 
based on the review of Grevesse & Noels (1993). The val- 
ues most relevant to our study are summarized in Col. 2 of 
Tab.[TJ A subset of models with solar-scaled abundances 
but log (Z/Zq ) = —0.3 was also computed but will only 
be discussed briefly in Sect. 16. 21 

The second set of models has the same metallicity 
Z=0.02 as the reference set, but modified abundances of 
4 He, 12 C, 14 N, 16 (Col. 5 of Tab.HJ. In the following, the 
adopted modified abundances will be refered to as "RSG- 
specific abundances" . 

The RSG-specific abundances were selected by inspec- 
tion of the evolutionary tracks of Schaller et al. (1992; their 
case of standard mass loss) for stars with initial masses 
above 7 Mq , at evolutionary timesteps with effective tem- 
peratures below 4500 K. The values selected are represen- 
tative of the final red supergiant stages of a star of initial 
mass 20 Mq (T off ~3550K). Stars of lower initial masses 
would have RSG abundances closer to the main sequence 
values, while the tracks at 25 M Q reach significantly larger 
modifications (tracks above 25 M Q don't extend to the low 
effective temperatures of red supergiants). Note that the 
initial mass fractions of Schaller et al. (1992) are not ex- 
actly the same as assumed in our reference set (mainly 
because of their larger He abundance), but that these dif- 
ferences are small compared to those that distinguish red 
supergiants from zero age main sequence stars. 

In Tab. [H the adopted abundances are compared to 
other values in the literature. The tracks of Meynet et al. 
(1994) assume larger mass loss rates than those of Schaller 
et al. (1992). More ingredients distinguish the models of 
Brcssan et al. (1993) from those of the Geneva group. 
Nevertheless, predicted surface abundance alterations are 
of a comparable amplitude. Evolutionary tracks for stars 
with initial rotation predict that comparable 14 N enhance- 
ments are reached already by the end of the main sequence 
(Meynet & Maeder 2000, 2003), and are further increased 



during the RSG phase. The achieved abundance ratios de- 
pend strongly on initial rotation velocity, on initial mass 
and on mass loss prescriptions. One major difference be- 
tween rotating and non-rotating models is the length of 
time a star spends as a red supergiant with modified 
abundances in one and the other case (see also Maeder 
& Meynet 2001). 

Mixing also occurs along the red giant branch and 
asymptotic giant branch for low and intermediate mass 
stars. The RSG-specific abundances adopted here are 
more extreme than those obtained through 1st dredge- 
up on the RGB (e.g. Iben 1964, Charbonnel et al. 1996, 
Girardi et al. 2000). In particular, the RSG-specific en- 
richment in He and 14 N and the drop in H and 12 C 
are larger than expected from 1st dredge-up. More mix- 
ing may however occur through second dredge-up on the 
early asymptotic giant branch (for stars with initial masses 
above about 4M Q , e.g. Becker & Iben 1979, Boothroyd 
& Sackmann 1999), and though "non-standard" extra 
mixing for low mass stars that have evolved on the 
red giant branch past the RGB-bump (Charbonnel 1994, 
Charbonnel & do Nascimento 1998). Both these processes 
affect relatively luminous giant stars. The second one 
seems to be less efficient at the quasi-solar metallicities 
we consider than for population II stars. Therefore, we 
expect our RSG-specific abundances not to be appropri- 
ate for most solar neighbourhood giants of class III, while 
they might be relevant to some giants of class II. 

We note that future calculations with modified surface 
abundances will include mixing-induced enhancements in 
the 13 C abundance, since 13 CO is a clear feature in near- 
IR spectra of cool stars. The effects of recent changes in 
measurements of the solar abundance ratios (Asplund et 
al. 2005) will also be investigated. 

2.3. Spectra in numerical form 

The model spectra for M=l Mq with solar abundances, for 
M=15Mq with solar abundances, and for M=15M Q with 
RSG-specific abundances, are made available in FITS for- 
mat through CDS. Because the quality assessments made 
in this paper are restricted to resolutions of order 10 3 
in the near-IR (and a few hundred at wavelengths below 
0.81 //m), the spectra made available are smoothed with 
a Gaussian to a full width at half maximum of 2 A. The 
initial models, calculated with a wavelength step of 0.1 A, 
can be requested from A.L. or P.H.H. 

3. Trends in the models 

Spectra illustrating the effects of mass, gravity and sur- 
face abundances are provided in Fig. [21 In this section, we 
will discuss quantitative trends using selected colours and 
molecular band indices. The indices are measured for each 
spectrum using: (i) the standard J, H, K filter passbands 
of Bessell & Brett (1988); (ii) narrow and intermediate- 
band filters as described in Tab. [5] All narrow and inter- 
mediate filter passbands are approximated with rectangles 
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Table 1. Surface abundances (by mass) 





Adopted 


Geneva 


Padova 


Adopted 


Geneva 


Padova 




reference 


1992-1994 


1993 


RSG-specific 


1994 


1993 


Element 


set 


ZAMS 


ZAMS 


set 


RSG 


RSG 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




L U 


0.703 


0.680 


0.700 


0.580 


0.55 


0.63 


4 He 


0.280 


0.300 


0.280 


0.400 


0.43 


0.35 


12 C 


0.0030 


0.0044 


0.0049 


0.0022 


0.0020 


0.0032 


14 N 


0.00092 


0.0014 


0.0012 


0.0065 


0.0080 


0.0051 


16Q 


0.0083 


0.0106 


0.0106 


0.0076 


0.0068 


0.0084 



Notes: Col. 2: Abundances adopted in our reference set of solar metallicity models. Col. 3: For comparison, main sequence 
abundances of Schaller et al. (1992), also used by Meynet et al. (1994). Metal abundance ratios based on Anders & Grevesse 
(1989). Col. 4: Main sequence abundances of Bressan et al. (1993). Metal abundance ratios based on Grevesse (1991). Col. 5: 
Abundances adopted in our RSG-specific set of models (see text). Col. 6: Final RSG abundances of Meynet et al. (1994) for 
20 M stars. Col. 7: Final RSG abundances of Bressan et al. (1993) for 20 M stars. 



Table 2. Filter and index definitions 



Filter 


Center 

(H 


Width 

(A) 


Notes 


104 


1.0395 


50 


quasi-continuum (1) 


108 


1.0800 


60 


CN (1) 


110 


1.1000 


50 


quasi-continuum near CN (1) 


220 


2.2000 


1100 


quasi-continuum (2) 


236 


2.3600 


800 


1st overtone CO (2) 


COH 


1.6222 


80 


2nd overtone CO 


COHcl 


1.6160 


30 


absorption-minimum near CO 


COHc2 


1.6285 


30 


absorption-minimum near CO 


Index 


Definition (3) 




104-220 
CO(2.3) 
CO(1.6) 
CN(l.l) 


-2.51og(104/220)+cst 

-2.51og(236/220)+cst 

-2.5 log [2 COH/(COHcl+COHc2)]+cst 

-2.51og(110/108)+cst 



Notes : (1) Adapted from the 8-colour system of R.F. Wing 
(White & Wing 1978). (2) Adapted from Frogel et al. (1978). 
(3) est stands for a constant that gives the index the value 
for a model spectrum of Vega. 



of identical central wavelength and width as the filters in 
the original references (as already done by Bessell et al. 
1989). A model Vega spectrum provides the zero points in 
all passbands. 

3.1. Colours 

As shown in Fig. [31 colours that combine flux measure- 
ments around 1.04 /im, in the J band and in the K band 
are good indicators of T eS in theory, as their sensitivity to 
surface gravity is low. Above 3400 K, the spread in log(g) 
corresponds to a full spread in T cS of about 200 K for the 
15 Mq models (left panel). For 1M models, the corre- 
sponding spread is much smaller: about 60 K, centered on 
a line very close to the models at 15 M Q and log(g)=0. 
At the lowest temperatures, contamination of the pseudo- 
continuum in the K band with H2O absorption leads to 
reduced fluxes in low gravity stars. Unfortunately, in the 
two-colour plots useful for observers the extinction vectors 
run almost exactly parallel to the temperature sequence 



(right panel) : more resolved spectral information is nec- 
essary to estimate an effective temperature from near-IR 
data. 

Figured] illustrates the gravity dependence of colours 
involving H band fluxes. At high gravities, the minimum 
of the opacity of H - around 1.6 /im produces a distinct 
hump in the H band spectra, with correspondingly blue 
H-K and red J-H colours. At low gravities, molecular ab- 
sorption due mainly to CO and CN erases this continuum 
opacity feature. Such an effect was already mentioned by 
Bessell et al. (1991), though their interpretation probably 
underestimated the role of CN as compared to CO. The 
observations of Lancon et al. (2007) and those described 
in Kucinskas et al. (2005) provide a convincing validation 
of the JHK colours of the new models. 

The effect of mass on the H band flux is insignificant 
at log(g)>0. For lower gravities, H-K increases by up to 
only 0.02 magnitudes when going from 15 to 1M , at 
T cff >4000K. 

Switching from solar-scaled to RSG-specific abun- 
dances has the following (small) effects on the above 
colours. All colours tend to become bluer. Colour dif- 
ferences in H-K and J-H remain smaller than 0.04 mag 
(and are <0.02mag for most stellar parameters). The 
colour index 104—220 (defined in Tab.[5]) is reduced by up 
to 0.08 mag. The bolometric corrections to the K band, 
BC(K), are essentially unchanged (the PHOENIX values 
agree with those of Levesque et al. 2005 to within a 
few percent between 3400 and 4300 K). Effects this small 
would be difficult to exploit on real stellar spectra. 

3.2. Molecular indices 
3.2.1. CO 

CO is a long known indicator of luminosity class (Baldwin 
et al. 1973). It is sensitive to gravity and effective temper- 
ature, but also to metallicity and micro-turbulence. As 
indicated previously, a constant micro-turbulent velocity 
of 2km.s _1 is used in this paper except in a few models. 
Large micro-turbulent velocities deepen the 1st overtone 
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Fig. 2. Effects of gravity, temperature and surface abundances on PHOENIX model spectra. In each figure, the upper 
spectrum has T off = 4000 K and the lower one T eS — 3300 K. For easier comparison, fluxes have been normalized to 
values comparable to those in the upper left diagram. The figures on the left are at log(g)=2, those on the right at 
log(g)=-0.5. The upper figures are for solar abundances, the lower ones for RSG-specific abundances. The effect of 
mass (1M© vs. 15 M©) is too small to be easily identified on this type of figure. 



band of CO more than the 2nd overtone band, because surface abundances are shown. As commonly found, CO 
line saturation is more important in the former than in increases with decreasing temperatures and gravities. The 
the latter (e.g. Origlia et al. 1993, 1997 and refs. therein). CO strength progressively levels off when log(g) takes neg- 
ative values (i.e. the further dependence on g is negligible). 
On the left panel of Fig.® the changes of the 1st over- Contamination by H 2 at high log( 5 ) produces a drop of 
tone CO band at 2.29 ^m with gravity, temperature and 
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Fig. 3. Temperature sensitive near-IR colours (solar abundances, 15 M Q ). In the right panel, the effect of extinction 
on cool stellar spectra is shown for Ay=l (using the extinction law of Cardelli et al. 1989, with Ry=3.1). 
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Fig. 4. Gravity sensitivity of H-K (solar abundances, 15 M©). Extinction vector as in Fig.[3j 




Fig. 5. Measurements of the strength of the 1st overtone CO band at 2.29 /im and of the 2nd overtone CO band at 
1.62 /im in the model spectra (15 M ). Symbols: temperature sequences at the indicated gravities, for RSG-specific 
abundances. Lines: corresponding sequences for solar-scaled abundances (black and light-coloured lines alternate, and 
the dashed line has log(g)=0). Note that at a given CO (2.3), xCO (1.6) tends to be weaker in low gravity stars. 
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the CO index below 3200 K. Switching from solar to RSG- 
specific abundances reduces the CO strength generally by 
small amounts. The effect is maximum around 4500 K in 
low gravity models : RSG-specific models with log((/)=— 1 
at 4500 K have the same CO index as solar abundance 
models with log(g)=2, or alternatively with log((/)=— 1 
but 4800 K. 

The effects of log(g), T eB and abundances on the ap- 
parent strength of the 2nd overtone CO band at 1.62 /im 
are similar, with two notable exceptions. First, the effects 
of changes in the abundance ratios are smaller than at 
2.29 /im. Second, low gravity saturation is reached earlier. 
The result is summarized in the right panel of Fig.[5j In 
particular, low gravity stars tend to have weaker CO bands 
around 1.6 fim than high gravity ones, at a given strength 
of the 2.29 [im band. Contamination of the H-band fluxes 
by CN absorption contributes to producing this trend, as 
already hinted at by Wing & Spinrad 1970. 

Moving down from 15 M Q to 1 Mq has a negligible 
effect on the near-IR CO bands for log(g)>0. At lower 
gravities, the 1M Q CO bands are weaker than the 15 M Q 
bands (the effect is stronger for the 2.3 /xm band than for 
the 1.6 /xm bands). 

3.2.2. CN 

CN displays prominent near-IR absorption bands that 
have been studied extensively in the context of carbon 
star models (e.g. Loidl et al. 2001). The CN bands are 
prominent in red supergiants as well (White & Wing, 
1978). While in carbon stars essentially all surface oxy- 
gen is locked into CO, CN coexists with other oxides in 
the atmospheres of red supergiants. 

The behaviour of CN bands with varying model pa- 
rameters is complex, as shown in Fig. [6] Bessell et al. 
(1989) describe the decrease of the CN 1.1 /xm band 
strength with decreasing effective temperature below 
3800 K, as well as its gravity dependence (stronger CN for 
lower gravities). Our more extended temperature range 
shows that the maximum CN strength is reached between 
4200 and 4800 K. Both the location of the maximum and 
its actual strength depend on surface gravity, and on the 
chemical composition of the atmosphere. CN bands are 
strongly enhanced in models with RSG-specific abundance 
ratios. The effect of mass is small. CN is also enhanced 
when larger micro-turbulent velocities are assumed (Tsuji 
1976). In empirical samples, spectra with strong CN ab- 
sorption bands compared to their CO bands are candi- 
dates for modified surface abundances. 

3.2.3. Other molecular bands longwards of 1/im. 

H2O. H2O appears very abruptly in the models below a 
gravity-dependent threshold temperature: 3600 K at 
log(<?)=l, 3100 K at log(<?)=— 1 (based on measure- 
ments in the K-band wings of the H2O band centered 
near 1.9 /im). Below this threshold T cff , higher gravi- 



ties lead to stronger spectral bands. Varying the mass 
between 1 and 15 M , or switching to RSG-specific 
abundances, has only small effects on the H2O bands. 

TiO. Near-IR TiO bands around 1 /im (TiO 5, Av = -1) 
and 1.25 /im (TiO <fi, Av = —1) appear progressively 
in the models below a gravity-dependent temperature: 
-3600 K at log(g)=l, -3400 K at log( 5 )=-2 (based on 
visual inspection of the spectra in that region) . Other 
near-IR TiO bands longwards of 1 /im, such as the 
0, Av — band near 1.12 /im are hidden in CN absorp- 
tion. Again, varying the mass between 1 and 15M Q , or 
switching to RSG-specific abundances, has only small 
effects. We note that the next version of PHOENIX cal- 
culations will include an update of the TiO partition 
function and of the electron /-values of the TiO bands 
for the AMES TiO line list (Schwenke 1998), which 
appears to improve spectral synthesis results for M 
dwarfs (Allard et al., in preparation). 

VO. The 1.05 /tm VO band (VO A-X Av = 0) is sig- 
nificant in the model spectra only at T eff <3200K for 
log(<?)<l- The effect of mass or abundances is small. 

4. Models versus data in two-colour plots 

In this and the following sections, we compare the models 
with the data collected by Langon & Wood (2000) and 
Langon et al. (2007). Both sets provide spectra at a spec- 
tral resolving power of order 1000 between 0.97 /im (some- 
times 0.81 /tm) and 2.4 /tm. The first set adds low resolu- 
tion extensions through the optical range down to 5100 A 
for a few of the stars. The merged sample contains lumi- 
nous stars of spectral types G3I to M5I, G3II to M3II, 
and G5III to M5III, as well as asymptotic giant branch 
variables for comparison. 

As shown in Fig.0 the agreement between models and 
data is excellent in near-IR two-colour plots once extinc- 
tion has been accounted for (see Sect. [5]). Note that in this 
and the following figures the lines join models at constant 
log(</), while the data should follow real red giant or red 
supergiant branches. Evolutionary tracks predict that the 
warmer giants have log(gr)>2 while the cool ones reach 
log((/)~0. Red supergiants of various luminosity classes 
are expected to have log(g)<0.5. 

Figure[8] combines measurements of the first and sec- 
ond overtone CO bands and the 1.1 /im CN band with 
J-K. Agreement between solar metallicity models and em- 
pirical data is good for giant stars. The strongest offset is 
noted for the second overtone CO bands in the H window, 
which tend to be too strong in the cool giant star models. 
The figures also suggest that modeled first overtone CO 
bands might be slightly too weak at low T cS . Because ex- 
tinction affects J-K, two- index diagrams with a negligible 
sensitivity to reddening are presented in Fig. [9] The same 
conclusions hold. The CO line list data are from Goorvitch 
& Chackerian (1994 a, b), and are known to work very well 
in the case of M dwarfs. Therefore, it is unlikely that the 
line list data for CO is the cause of the CO band discrep- 
ancies. We note that 13 CO contributes to the measured 
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Fig. 6. Measurements of the strength of the CN band at 1.1 /im in the model spectra (15 M ). Symbols and lines are 
as in Fig.O The main effect of a decrease of mass (15— >1 M Q ) is a shift of the models at log(<?)<0 to the left by up to 
0.04 mag in the right hand panel. 
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Fig. 7. Two-colour plots with observational and calculated data. The thin lines are T eS sequences for solar abundances, 
the thick lines are for RSG-specific abundances. Solid lines are at log(<?)=2, dashed lines at log(g)=— 1. The dotted line 
follows models along an illustrative red giant branch at solar metallicity ([T cff (K),log(g)] = [4200,2], [3800,1], [3400,0]). 
Symbols as explained are measurements on dereddened versions of the spectra of Langon et al. 2007 (see Sect. [5]), and 
dots O-rich Miras from Langon & Wood 2000 (not dereddened). The reddening vectors are as in Fig. [3] 



strength of the first overtone CO bands, and that changes 
in the 13 C abundances induced by stellar evolution may 
be responsible for some systematic effects. Slightly larger 
micro-turbulent velocities could improve the ratio of the 
first to the second overtone CO band strengths, but would 
also affect the CN bands. The outlier giant star near J- 
K=0.7 with weak band strengths is the only Population II 
star of the observed sample (HD 218732). Eye inspection 
of its spectrum reveals a metal poor atmosphere immedi- 
ately. By contrast, the other giant stars appear to form a 
reasonably homogeneous sample in terms of metallicity. 

While the trends with gravity present in the predicted 
molecular bands agree qualitatively with those observed, 
the molecular bands of only a fraction of the observed red 
supergiants can be reproduced quantitatively. Models with 



RSG-specific abundances are favoured for a significant 
number of these objects, which show stronger CN bands 
than the solar metallicity models can produce. However, 
the CN and CO measurements show that despite the im- 
provement achieved with the adopted changes in surface 
abundances, the model parameters explored here are not 
able to account for the whole range of molecular band 
strengths observed. Models with larger micro-turbulent 
velocities reach further into some of the areas occupied by 
real red supergiants, justifying the ongoing extension of 
the model grid. Alternatively, some red supergiants in na- 
ture may require even higher 14 N abundances than tested 
here or effective gravities lower than log(g)=— 1. 

Comparisons between models and data on the star-by- 
star basis are provided in the following section. 
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Fig. 9. Molecular index plots with dereddened data. Symbols and lines are as in Fig. [7] 



5. Direct comparison between observed and 
theoretical spectra 

5.1. Method 

The comparison between models and data is based on 
the computation of reduced \ 2 differences. The theoretical 
spectra are smoothed to the resolution of the data, using 
gaussian broadening functions with adequate widths (note 
that for the optical spectra, whose resolution was seeing- 
dependent, we adopt a single smoothing value which in 
some cases underestimates the actual AA). They are then 
resampled to the wavelength step of the data, i.e. 5 A. A 
window function is used to eliminate the spectral intervals 
most strongly affected by telluric absorption, around 1.15, 
1.4 and 1.9 fmx. 

The rms noise of the data is modelled as being propor- 
tional to the square root of the signal. Numerical values 
given below assume an average signal-to-noise ratio of 50. 
This simple noise model is a reasonable representation of 
the typical high frequency noise of the data. The addi- 
tional uncertainties due to flux calibration errors are not 
explicitly accounted for. They lead mainly to uncertain- 
ties in the estimated extinction values. A further discus- 
sion of the effects of the weighting of the data is provided 
in Sect.EU 

A mass of 1 M Q is assumed for giants and bright giants, 
a mass of 15 M Q for supergiants (see Sect.[|5]). For each 
empirical spectrum, the adopted algorithm loops through 



all model temperatures and gravities (separately for the 
two sets of abundances). At each of these points, it also 
loops through an adequate set of extinctions (using the 
extinction law of Cardelli et al. 1989 with Ry — 3.1), and 
minimizes the \ 2 with respect to the extinction parameter 
Ay. The step in Ay is of 0.05 mag. A \ 2 map is produced 
in T eB —log(g) space, and the 9 best fits (of the T eB — log(p)- 
Ay space) are plotted for inspection. The x 2 value of the 
9th best fit is typically higher than for the 1st best fit 
by 10 % when the best fits are good, and by only a few 
percent when they are poor. Typical uncertainties on the 
derived parameters in the case of good fits are ±100K 
in T cB , and ±0.2 in Ay. For these good fits, gravity is 
usually determined to better than one step within our set 
of models (log(<?)=— 1, — 0.5, 0, 1, 2). Preliminary models 
with micro-turbulent velocities larger than 2kms~ 1 were 
tested only for supergiant star spectra and, in this paper, 
they are only discussed in cases where the initial fits were 
poor. 

The method is robust with respect to errors in the po- 
sitions of the individual molecular lines that jointly define 
the pseudo-continuum of the spectra at the resolution of 
interest here. Such errors were noted as a difficulty in the 
measurement of individual metal lines by Vanhollcbckc et 
al. 2006. In order to verify this, we added random noise 
to the data at the level of a few percent, i.e. enough to 
completely alter the apparent positions of the blended CN 
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Fig. 10. Fit to a late type giant star spectrum 
(HD 145480, M5III). The data are shown as dotted lines, 
the model as a solid line. The window function for the x 2 
fit is also shown. Model parameters are 3400 K, log(g)=0, 
Ay=0.4, x 2 = 1-7. Such a fit quality is typical for all the 
giants without available optical spectra. 



lines between 2 and 2.3 fiia. Differences in the derived pa- 
rameters were well within the uncertainties stated above. 



5.2. Giant stars 

The data sample contains 35 stars of class III with spec- 
tral types G5 to M5 (after elimination of one particularly 
uncertain luminosity class and one known metal poor star 
with obviously weaker spectral features). Their shortest 
wavelength is 0.97 /im, except for 5 spectra with an opti- 
cal extension. 

Good fits are obtained for essentially all the near-IR 
spectra with the solar metallicity models. An example is 
given in Fig.[TUJ Among the satisfactory fits, there is a 
tendency for \ 2 to increase systematically from about 1 
for types < K4 (x 2 =0.6-1.4, depending on the actual S/N 
of individual spectra) to about 2 (1.5-2.5) for the M stars. 
This trend is due to a wealth of weak lines and deeper 
molecular bandheads at low temperatures, which among 
others induces a higher sensitivity of the value of the x 2 to 
residual wavelength calibration errors or model line lists. 

The x 2 values for combined optical+near-IR spectra 
take values of 3 to 7 for satisfactory fits (considering the 
sensitivity of the x 2 to the smoothing parameter at op- 
tical wavelengths and to flux calibration errors over such 
an extended wavelength range). Examples are shown in 
Fig.[TTJ The best fit to the spectrum of the coolest giant, 
the M5III star BS 4267 (= HD 94705) requires marginally 
negative extinction. Considering the flux calibration un- 
certainties and the choice of one particular extinction law, 
such a result is not alarming. Cooler models by only 100 K 
or models with higher log(g) by 0.5 would result in a pos- 
itive value of the estimated Ay. The most obvious short- 
coming of the models for giants this cool is an overpredic- 
tion of the TiO absorption bands near 1 and 1.25 /im. 




1.5 

Wavelength (/xm) 

Fig. 11. Best fits to four giant star spectra that extend to 
optical wavelengths. Data are shown as dotted lines, best 
fit models as solid lines. The window function used to re- 
ject the noisiest spectral regions is also shown. From top to 
bottom: K3.5III star BS 4432 with model [T eS , log(g), Ay, 
X 2 ] = [4100 K, 2, 0.55 mag, 2.3]; M0III star BS4371 with 
[3900 K, 1, 0.2 mag, 4.5]; M2III star BS 5301 with [3800 K, 
1, 0.55mag, 5.2]; M5III star BS4267 with [3300K, 1, - 
0.3 mag, 13.9]. 



The results of the fitting procedure can be summarized 
as follows: 

Temperatures range from 5300 K for type G5, to 
3300 K for type M5. As expected from stellar evolution 
tracks, the highest available gravity (log(g r )=2) is selected 
for giants earlier than K7 (with one exception), then pro- 
gressively more spectra are assigned log(<?)=l and later 0. 
Ay values are spread between and 1 (with 4 cases of 
marginally negative values). No correlation is found be- 
tween Ay and T eff . 

Adopting the models with RSG-specific abundances 
rather than solar ones leads to poorer fits in all but 
one case. The values of the reduced \ 2 increase by 0.5 
units on average. While the distribution of estimated 
effective temperatures for the sample is relatively flat 
with the assumption of solar abundances, it becomes bi- 
modal with RSG-specific abundances : temperatures be- 
tween 4500 and 5000 K are systematically avoided, be- 
cause the CN bands of these models are too strong at 
the surface gravities of giant stars (cf. Fig. [6]). This result 
was expected from Figs. and [H It is satisfactory as our 
set of RSG-specific abundances is not designed to match 
abundances in giant stars. 
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Fig. 12. Best fits to four bright giant star spectra (class 
II) that extend down to 0.81 /an. Figure set-up is as in 
Fig.fnj From top to bottom: G8II star HD 150416 with 
[T eff , log( ff ), Ay, x 2 ] = [5000K, 1, 0.1 mag, 0.65]; K2II 
star BD-29 2374 with [4500 K, 2, -0.15 mag, 0.65]; K5II 
star HD 168815 with [4100K, 1, 1.3mag, 1.8]; M0.5IIb 
star HD 132933 with [4000K, 2, 0.4mag, 1.4]. Among the 
above, we classify the K5II fit as satisfactory, the others 
as very good. 
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Fig. 13. Best fits to two bright giant star spectra for 
which the CN band at 1.1 fim is particularly poorly re- 
produced, and of the coolest class II star of the sample. 
Figure set-up is as in Fig. QT] From top to bottom: M0Ib/II 
star HD 145384 with [T cff , log( 5 ), Ay, X 2 ] = [3400K, - 
0.5, -0.15 mag, 3.2]; M0.5II star HD 142676 with [3900 K, 
2, 0.0 mag, 4.0]; M3II star HD 153961 with [3500 K, 0, 
1.0 mag, 3.3]. 



5.3. Bright giants 

The sample contains 29 bright giants of class Ib/II or II. 
Spectral types range from G3 to M3. None of the corre- 
sponding spectra extend through optical wavelengths, but 
11 extend down to 8100A. Their properties in terms of 
colours and molecular indices are spread between those of 
giants and supergiants. On average, they display slightly 
stronger bands of CO and significantly stronger bands of 
CN than giants of class III, at a given (dereddened) colour. 

The solar metallicity model fits to all the spectra are 
satisfactory, two thirds of them being very good (Fig. [T2|) . 
The models clearly contain all the molecular bands re- 
quired. Marginally negative values of Ay are obtained in 
four cases, which is again not unexpected considering flux 
calibration uncertainties. The most common shortcomings 
found when the fits are not perfect are the following: 

- There is a tendency for the models to show stronger 
CO bands at 2.29 /im and weaker CN bands at 0.93 and 
1.1 fim than observed. This problem is only detectable 
clearly when the data extend down to 0.81 /im. 

- For stars with spectral types around K5II whose spec- 
tra extend down to 0.81 /im, the models struggle to repro- 
duce the energy distribution around 1 /im, where it peaks 
between deep CN bands. This difficulty is certainly related 
to the strength of the CN bands at those temperatures (see 
Fig.EJ. 

- In two cases (spectral type M0Ib/II and M0. 511), the 
model CN bands are too weak while CO is reproduced 
well (Fig.HU). 

The T cS and log(<7) scales obtained for the bright 
giants have a scatter similar to those found for giants 
and supergiants and are located between the two, as ex- 
pected. Bright giants with spectral types earlier than K4 
(included) are assigned log(<?)=2 (with one exception: 
HD 170457, G8Ib/II), and values of \og(g) for types K5- 
M3 are scattered between 1 and -0.5. No correlation is 
found between Ay and T cS . 

When moving from solar abundances to RSG-specific 
abundances, the x 2 test indicates that the fits are de- 
graded in a majority of cases (typically by 0.5 \ 2 units, as 
for the sample of class III stars). However, a significantly 
improved x 2 is obtained for 4 stars, and the x 2 changes 
are insignificant in 7 cases. 

The improvements correspond to four stars of type K2 to 
M0 (out of the 11 bright giants available in this range), 
with estimated T eS of 4300 to 3400 K. Eye inspection of 
the corresponding four spectra shows that the decrease in 
X 2 corresponds to a better fit to the CN bands, which were 
not deep enough (by small amounts) in the solar metallic- 
ity models. In some cases, the improved x 2 was associated 
with a decrease in T eff by 100 K or an increase of log(g) 
by one bin size, but statistics are too small to define sig- 
nificant trends. 

Degraded fits are frequently associated with excessive 
strengths of the model CN bands when the RSG-specific 
abundances are used. The T eff -distribution obtained with 
the RSG-specific abundances still shows a zone of avoid- 
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Fig. 14. Good and reasonable best-fits to warm red su- 
pergiant spectra (class I). Figure set-up is as in Fig.fTTI 
Abundances are solar. From top to bottom: G2Ib star 
HD 182296 with [T cff , log( 5 ), Ay, x 2 \ = [5000K, 0, 
0.95 mag, 1.3]; G8Iab star HD 206834 with [4900 K, 1, 
0.3 mag, 0.85]; G5Iab star HD 170234 with [4500 K, 0, 
2.15mag, 1.8]; K3Iab star HD 187238 with [4100K, 0, 
1.65 mag, 2.2]. We have not counted the K3Iab case as 
a good fit, because the model CO bands around 1.7 /jm 
are too strong. Note that the spectral type of the second 
and third stars are likely to be incorrect. 



ance between 4500 and 5000 K, but the effect is not as 
obvious as in the case of class III stars. Although small 
number statistics affect this result, we note that all class 
II spectra with estimated T off in that range have poorer fits 
with the RSG-specific abundances than with solar ones. 
As expected, models with the adopted RSG-specific abun- 
dances do not apply to the majority of class II stars but 
they do allow us to identify candidate objects that may 
have suffered more than standard dredge-up. 

5.4. Supergiant stars 

The data sample contains 37 spectra of stars of class I, la, 
lab or lb (after removal of one particularly odd case that 
is probably misclassified and of one spectrum with poor 
correction for telluric absorption). Spectral types range 
from G2 to M5. The stars with the latest spectral types 
(> M2) are all known or suspected variables (as are the 
vast majority of late type supergiants in nature). 9 spec- 
tra, all with spectral type M, extend through the optical 
range; note that the optical and near-IR spectra of indi- 
vidual objects were taken within less than three weeks of 
each other. 8 spectra extend down to 0.81 //m. 
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Fig. 15. A selection of marginally acceptable best-fits to 
red supergiant spectra. Figure set-up is as in Fig. [141 
From top to bottom: G5Iab star HD 165782 with [T eS , 
log(g), Ay, x 2 } = [4900 K, -1, 2.4 mag, 2.8]; MHab star 
HD 98817 with [3700K, -1, 1.8 mag, 7.8]; M2Iab star 
BS3364 (=HD 72268) with [3500 K, -0.5, 1.0 mag, 12.4]; 
M4.5I star V774Sgr with [3200K, -0.5, 1.7 mag, 22.7]. 



Good fits to the red supergiant spectra with solar 
metallicity models are obtained for 13 of the 37 spectra, all 
of which are of spectral type G2-G8. 16 of the remaining 
spectra find a model representation that is still reasonable 
(though often significantly poorer than the fits we called 
satisfactory within class II above). These are spread over 
the whole range of spectral types and include some of the 
data that extend through optical wavelengths. In general, 
stars of luminosity class lb are easier to fit than those of 
class I, la or lab, and all class lb stars of our sample can 
be matched well or reasonably well. Finally, we classify 7 
of the red supergiant fits as poor. Five of these correspond 
to variable stars with spectral types later than M3.5 (class 
I, la or lab), the two others are of spectral type MOIa and 
G5Ia. 

Figures Q4] [15] and [16] illustrate some of the good, in- 
termediate and poor model fits. 

The main shortcomings of the models when the fits are 
of intermediate quality are the following : 

— A relatively common feature is a shortage of flux in the 
models around 1 /nm, as seen in two spectra of Fig.Qj)] This 
problem was already mentioned for a few bright giants of 
class II, as a property that is associated with strong CN 
bands and can be identified only when the observed spec- 
tra extend far enough to short wavelengths. 

- Even when the 1st overtone CO bands after 2.29 /im 
are reproduced reasonably well, it happens that the rel- 
ative strengths of the 2nd overtone CO bands in the H 
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Fig. 16. A selection of poor best-fits to red supergiant 
spectra. Figure set-up is as in Fig.[TJ] From top to bot- 
tom: G5Ia star HD 155603 (classified KOO-Ia by Keenan 
and McNeil, 1989) with [T eB , log(g), Ay, X 2 ] = [4300 K, 
-0.5, 1.8 mag, 8.6]; MOIa star Trumpler 1-27 with [4200 K, 
-1, 4.4mag, 13.4]; M3.5I star IRC -20427 with [4000K, 
-1, 5.35 mag, 15.2]; M4-5Iab star CLCar with [3300K, 2, 
1.65 mag, 21.6]. The fits being of poor quality, the derived 
parameters are not reliable and are given here only for the 
sake of completeness. 

window are incorrect, the transitions at longer H window 
wavelengths (1.65-1.75 /im) being too strong in the models 
compared to the data (last spectrum of Fig. [14] and 2nd 
and 3rd spectrum of Fig.[T5]) . 

— In the coolest models, bands of TiO appear (in partic- 
ular near 1.25 /im) that are not seen in the data. 

Poor fits are obtained for the coolest stars (e.g. bot- 
tom spectrum of Fig. fTBj) and for stars with extreme CN 
bands (e.g. top three spectra of Fig. [16]). We recall that 
the coolest stars are also variable and that discrepancies 
are to be expected in a comparison with static models. 
When the CN bands are strong, the derived temperatures 
are a compromise between the necessity to reproduce the 
energy distributions and the CO bands at 2.29 /jm (which 
pulls towards low temperatures), and the need to maxi- 
mize CN depths (which pulls towards 4100 K, cf. Fig.[6]). 
When optical spectra are taken into account, the relative 
weight of the CN bands is reduced compared to CO, opti- 
cal features and the energy distribution. On the contrary, 
when only wavelengths between 0.97 and 2.4 /im are avail- 
able the role of the CN bands is large. This explains why 
in Fig.fTHlsuch a large difference in T eS is obtained between 
the M3.5Ia star (no optical spectrum, best fit T eff = 4000 K) 
and the M4-5I star (optical spectrum available, best fit 
T cff =3300K). The temperatures of the MOIa and M3.5I 



stars of that figure are most probably overestimated. For 
a similar reason, the temperature of the G5Ia star at the 
top of the figure may be underestimated (compare with 
the G5Iab star in Fig. [T5]) . 

A typical problem with the best fit models for the 
spectra with very strong CN is the relative strength of 
the various CO bands. The G5Ia star HD 155603 (Fig.QU) 
provides the most extreme example. It has the strongest 
2.29 /im CO band of our whole supergiant sample and 
among the strongest CN bands as well, but in the H win- 
dow CO is almost absent. None of the current models with 
Vmic = 2kms _1 reproduces this combination. Models with 
larger micro-turbulent velocities improve the representa- 
tion of these extreme spectra. 

Water bands are another cause of disagreement be- 
tween models and data. The near-IR bands of H2O and 
CN overlap in wavelength to such a degree that confu- 
sion can occur (and has occurred in the past, cf. Wing & 
Spinrad 1970). The shapes of the 1.1 /im bands of H2O 
and CN are subtly different (cf. Fig. 5 of Langon & Wood 
2000). The bands observed in red supergiants correspond 
closely to CN, although contamination with H2O is pos- 
sible. The H2O band around 1.9 /im, which is very deep 
and broad in Miras, is inconspicuous in red supergiants. 
It may be present at a low level in the coolest ones ob- 
served, such as CL Car (Fig. [T6"l) , which are semi-regular 
long period variables. The clearest H2O signature in the 
observed red supergiant spectra is a sharp bandhead at 
1.32 /im, although the detection of this feature requires 
good corrections for telluric absorption bands. Based on 
this signature, the 7 coolest supergiants of our sample con- 
tain H 2 (all these are variable). The models however ei- 
ther do not show this bandhead (low g) or, when they do 
(high g), also display a 1.9 /im band that is much wider 
and deeper than observed. 

Finally, the semi-regular variables V774Sgr, EV Car 
and CL Car (Figs. [TBI and [T6| have a clear VO absorption 
band at 1.05 /im and small or inexistent absorption bands 
at 1 /im and 1.25 /im, two properties that are not matched 
by the models. 

When moving from the models with solar abundances 
to the RSG-specific abundances, the \ 2 test indicates that 
about a third of the fits are improved, another third are 
degraded, and the quality of the final third of the fits is 
essentially unchanged. 

The deteriorations, when present, are not severe. In 
most cases, it seems that abundance values intermediate 
between the adopted solar and RSG-specific sets would 
provide optimal fits, which is not surprising considering 
that evolutionary tracks for red supergiants cover a range 
of abundances. Eye inspection shows that quite a few stars 
with equally good fits with both model sets also fall in this 
category. 

The improvements obtained with RSG-specific abun- 
dances for a fraction of the red supergiants are signifi- 
cant, although they clearly do not resolve all the difficul- 
ties. They are associated with a better representation of 
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Fig. 17. Top: K4Ib star HD 185622a and best fit model 
with RSG-specific abundances. Middle: Residuals of the fit 
shown above (data— model). Bottom: Residuals of the best 
fit with solar metallicity models (x 2 1-22 times larger than 
in the RSG-specific case). Note the CN residuals below 
1 /zm, around 1.1 /im and in the slope between 1.45 and 
1.75 fxui. These are typical and systematic shortcomings 
of the solar metallicity models in the cases where RSG- 
specific models provide a better fit. 



the observed CN bands and sometimes also with a better 
match to the CO bands around 1.6 /zm (see also Sect. l5.5j) . 
One may distinguish two subcategories of improvements. 
On one hand, some of the stars that already had rea- 
sonable model counterparts with solar abundances have 
better, often good adjustments with RSG-specific abun- 
dances. These are mainly stars of type G and K. An ex- 
ample is given in Fig. 1171 On the other hand, the improve- 
ments refer to stars that had poor fits with solar abun- 
dances, and for which the RSG-specific abundances lead 
to somewhat better but still unacceptable fits. These are 
the same 7 stars as mentioned earlier. The models can- 
not simultaneously reproduce their CO bands (1.6 and 
2.29 /im), their CN bands and their energy distribution. 
More extended model grids are needed to characterize 
these objects. Problems related to H2O, TiO and VO, 
when present, remain essentially unchanged. 

The explored changes in surface abundances induce 
changes in the best-fit parameters for the sample of ob- 
served stars with maximal amplitudes of ±200 K. For the 
sample as a whole, there is no strong correlation between 
the change in T cft and the actual value of T eff , which is not 
surprising considering that many fits are imperfect and 
that the behaviour expected from theory is complex (see 
Sect.EHl. 



The T off distribution of the red supergiant sample 
obtained under the assumption of RSG-specific abun- 
dances shows no anomaly. Scrutiny of the 2D distri- 
bution of estimated parameters in the log(<7)-T ofI plane 
suggests that a narrow zone extending diagonally from 
[log(s)=0,T cff =4000K] to [log( 5 ) = l,T cft =5000K] (with no 
extension to lower gravities) might nevertheless be under- 
populated. The statistical significance of this gap is low be- 
cause of small sample numbers. Its presence would favour 
a general picture in which RSG-spccific abundances are 
only relevant to red supergiants with large initial masses 
or to late stages of red supergiant evolution. 

5.5. Effects of the weighting of various parts of the 
spectra 

Because there are systematic differences between the best 
fit models and the observed spectra, the best fit model 
parameters depend on the weights given to the various 
spectral features in the fitting procedure. Our standard 
method weights the data based on a reasonable simpli- 
fied model for the high frequency noise in the data. This 
adopted weight is inversely proportional to the square root 
of the signal, i.e. spectral regions with large fluxes con- 
tribute more to the y 2 than regions with small fluxes. 
Since the spectra of cool stars peak around 1 /im (in the 
flux density units adopted in this paper), molecular bands 
near this wavelength are important in the fit. In prac- 
tice, this weighting makes CN bands relatively important 
and CO bands (around 1.6 /im and 2.3 /im) comparatively 
unimportant. 

If the noise was indeed gaussian, uncorrelated be- 
tween pixels, and exactly of the amplitude assumed, then 
our procedure would select the models with the largest 
likelihoods. This is not the case (flux calibration errors, 
wavelength-dependent gains and contributions of the read- 
out noise, etc.), and therefore our weighting is in some 
ways non-optimal. We may choose various alternative 
methods (see Decin et al. 2004 for an interesting discussion 
of comparable issues). First, we may decide to fit measure- 
ments of the depths of one or several features rather than 
spectral segments. Unfortunately, the selection of one or 
the other feature remains somewhat arbitrary. Second, we 
may decide to focus on either the optical or the near-IR 
spectral range. This circumvents the difficulty of reproduc- 
ing the global energy distribution (possible uncertainties 
in the relative flux calibration of the optical and near-IR 
data, uncertainties in the adopted extinction law, etc.). 
Third, we may keep the whole spectrum but explore the 
effects of alternative weightings. We briefly summarize be- 
low the main trends found while investigating these three 
options. 

Figures [18] and [19] show how three important near-IR 
molecular features in the models compare with their ob- 
served counterparts, when our standard weighting proce- 
dure is applied to select the best fit. Note that the corre- 
sponding figures for H-K, J-H, J-K or 104-220 (not shown) 
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Fig. 18. Top: Strength of the 2.3 fim CO band of the 
best fitting solar metallicity models versus strength of 
this band in the dereddened observed spectra (cf. Tab[2|). 
The dotted line highlights the one-to-one relation. Middle: 
Same figure for the adopted measure of the 1.6 /im CO 
band. Bottom: Same figure for the adopted measure of 
the 1.1 /im CN band. 



are very well behaved, which only states that the adopted 
extinction law is capable of dealing with the actual extinc- 
tion (and with flux calibration errors) rather well. 

As expected, systematic drifts away from the perfect 
match are smaller for the 1.6 fim CO features than for the 



Fig. 19. Same as the top and bottom panels of Fig.[T8lbut 

for supergiants only, and using models with RSG-specific 
abundances. 



2.3 fim CO band, the latter being located in a region of 
low flux (small contribution to the \ 2 )- The best- fit mod- 
els have systematically deeper 2.3 fim CO bands than the 
data for warm stars (types G and K), but systematically 
too shallow 2.3 /zm CO bands for the coolest stars (type 
M). By changing the weights in the fitting procedure (e.g. 
by assuming a constant signal-to-noise ratio), the 2.3 /im 
bands can be reproduced better, but at the cost of a loss 
of the fit quality at shorter near-IR wavelengths. 

The CN bands are reproduced well for giant stars. But 
they are too shallow in the best fit models for some of the 
bright giants and for the supergiants. Here, changing the 
fitting weights has a small effect compared to more fun- 
damental model parameters such as abundances, gravities 
or micro-turbulence. RSG-specific abundances move the 
bulk of the red supergiants into a satisfactory location 
(Fig.[[n]). With RSG-specific abundances, the fits to CO 
bands around 1.6 /^m are not fundamentally improved or 
degraded on average, while the first overtone CO bands 
(2.3 /im) of the best fits become shallower, i.e. too shal- 
low. By assigning CO more weight in the fits, it is possible 
to reduce this discrepancy while still observing the global 
improvement for CN. But with the current grid of models, 
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no fully satisfactory solution can be found for any weight- 
ing scheme. 

The weights given to various spectral ranges impact 
on the estimated stellar parameters. Examples have al- 
ready been given in Sect. 15.41 and further discussions can 
be found below. 

6. Discussion 

Providing estimates of fundamental stellar parameters is 
a major application of theoretical spectra. Our discussion 
focuses on the determination of T cS from near-IR spectra 
using the new PHOENIX models. 

6.1. Stellar mass 

We have mentioned in Sect. [3] that the effects of mass on 
colours and molecular indices, at a given T eS and log(<?), 
are small. The comparison between the best fit parameters 
obtained assuming 1M© and 15 M nevertheless reveals 
a trend worth highlighting : for stars with high surface 
gravities (log(g)=2), the temperatures obtained with 1 Mq 
models are systematically lower by ^100 K than those ob- 
tained with 15 Mq models. This is particularly relevant to 
giants of class III, but also to the warmer of the bright 
giants of class II. Unfortunately, we found no systematic 
differences between the x 2 values obtained with one or the 
other assumption on mass. Thus, it is not currently pos- 
sible to determine mass using spectral fits such as those 
performed in this paper. Mass has to be fixed a priori by 
other means. 

For luminous giants and supergiants, i.e. stars with low 
gravities, we found no systematic effect of mass on best- 
fit T cff or log(<7). The differences in T eS between the two 
assumptions are scattered around with typical values 
of ±100 K (more in cases where even the best fits are not 
satisfactory). We note a correlation between the difference 
in T eS and the difference in log(<?) : when a change in the 
assumed mass leads to a rise in the best-fit T cff , it generally 
also produces a rise of the best-fit value of log(g). 

6.2. T cH -systematics related to surface abundances: 
model predictions 

The effects of surface abundance ratios on T cfi estimates 
(and on derived gravities) are of larger amplitude than 
those of mass, and we therefore describe them in more 
detail. They can be studied by fitting a sample of solar 
metallicity models with models with RSG-specific abun- 
dances, using the procedure described for fits to observa- 
tional data. The results depend on the wavelength range 
adopted in the analysis. They are illustrated in Fig. [20] 
The amplitude of the effect is of several hundred Kelvin. 
If we call <5T cff the difference between the input and out- 
put temperatures (output minus input), we find no simple 
linear correlation between 8T cS and T eS . 
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Fig. 20. Effects of the surface abundances on estimates 
of T e ff and log(g). The input T e g and log(<?) refer to 
solar metallicity models. The output parameters to the 
values obtained when fitting the solar metallicity spec- 
tra with models with RSG-specific abundances, using the 
procedure described in Sect. 15. II Solid: Input log(<?)=— 1. 
Dashed: Input log(<7)=l. 



The figure based on near-IR data (0.97-2.4/im, with 
the window and weight functions of Sect. 15. lj ) is tightly re- 
lated to the behaviour of the near-IR CN bands. Output 
T eff values around 4400-4900 K (depending on gravity) are 
avoided because the CN bands of those RSG-specific mod- 
els are too strong. A similar effect is present when optical 
wavelengths are used (0.51-0.97^m), but it is combined in 
the low- T eff range with a variety of effects due to oxides. 
The difference between the optical and near-IR tempera- 
tures is largest between 3500 and 4200 K, where the fluxes 
below 0.75 /im transit rapidly from being almost nil to be- 
ing large. Eye inspection of the fits shows that the best fit 
models sometimes deviate wildly from the "data" in the 
range not included in the fitting procedure, while over the 
range really used fits are dangerously good. When optical 
and near-IR spectra are used jointly, compensations occur 
and the correct Teff is recovered (to within ±100 K) be- 
low 4200 K. Positive offsets of up to 400 K however persist 
above this temperature for all gravities. 

The output log(g) equals the input log(p) below about 
4300 K when using near-IR data, and is higher by one 
log(g)-sampling step at higher temperatures. When using 
optical data, the behaviour depends more strongly on the 
actual value of the input log(<7). For high gravities, 8\og{g) 
is positive (one sampling step) at T eS > 3600 if. For low 
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Fig. 21. Effective temperatures derived from fits to near- 
IR spectra (A > 0.97 ^m), compared with trends in the 
literature. RSG-specific abundances are used for class 
I stars, solar abundances for classes II and III. Solid 
lines: temperature scale for giants, from van Belle et al. 
(1999) for T off <5000K and from Schmidt-Kaler (1982) 
for T off >5000K. Dotted line: temperature scale for su- 
pergiant stars from Levesque et al. (2005). Default spec- 
tral types from the SIMBAD database (operated at CDS, 
Strasbourg, France) are used for this figure. 



gravities, 5\og(g) is nil at the lowest and highest temper- 
atures, but peaks with a value +2 around 4200 K. 

We note that corresponding plots can be produced 
for the "extinction" -correction Ay (which accounts for 
colour changes in the analysed wavelength ranges reason- 
ably well). The qualitative aspects of the graphs for 5 Ay 
are similar to those of 5T e g, with a maximal amplitude of 
±0.6 magnitudes. 

For comparison, we have performed a limited explo- 
ration of the effects of metallicity (with solar scaled abun- 
dance ratios) on the derived temperatures. Models at 
log(Z/Z Q ) = —0.3 were computed for log(g)=l and —0.5, 
and best fits to these were obtained using solar metallicity 
models. Plots similar to those in Fig. l20l were constructed. 
The effects of the change in Z on the derived T off is no- 
tably smaller than those just described for modified abun- 
dance ratios. When using optical wavelengths, the trend 
expected from the well known metallicity-temperature de- 
generacy is found (lower temperatures are required at 
lower metallicity to produce similar optical band depths). 
The offset varies between 100 K (low T cff ) and 200 K (high 
T cft ). At near-IR wavelengths, the correct temperatures 
are recovered unchanged except for a few deviations of 
±100 K. In both wavelength ranges, however, gravities 
higher than input are derived (by one gravity bin). For 
complementary discussions on metallicity effects, we refer 
to Kucinskas et al. (2006) and Levesque et al. (2006). 



Fig. 22. Same as Fig.[2lJ but using only spectral types 
from Keenan & McNeil (1989) or Buscombe (1998). Using 
solar abundances moves the G6 supergiant in this figure 
down 200 K, the K0 supergiant up 100 K, and the M su- 
pergiants above 3600 K up by 100 to 200 K. 



6.3. T cff estimates for real stars 

Figure[2T] compares the effective temperatures derived in 
this paper from near-IR spectra with temperature scales 
from the literature. For giants, the plotted reference scale 
(below 5000 K) is based on angular diameter measure- 
ments (van Belle et al. 1999). The number of red super- 
giants with angular diameter measurements is small. For 
supergiants, we therefore show the scale recently obtained 
from fits of MARCS model atmosphere spectra to optical 
spectra by Levesque et al. (2005). The agreement is good, 
but the scatter is large. 

As a sanity check, we may restrict our data sample 
to stars with optical spectra and discard near-IR wave- 
lengths (A > 1 /im) in the fitting procedure. This pro- 
vides temperatures that are a priori more directly re- 
lated to spectral types. In addition, we keep only stars 
with MK spectral types from Keenan & McNeil (1989) or 
Buscombe (1998), and with small variability (according to 
the Simbad database information). Using solar abundance 
models for direct comparison with the results of Levesque 
et al. (2005), we find that 8 of the 9 stars in the subsample 
have estimated temperatures within less than 50 K of the 
reference relations. Most of the stars in the subsample are 
supergiants and all are of type K5 or later. Thus, in this 
range of parameters, there is no indication of a systematic 
difference between the temperatures derived from optical 
spectra using the new PHOENIX models or the MARCS mod- 
els of Levesque et al. (2005). 

To illustrate what fraction of the scatter in Fig.l2"T1mav 
be due to spectral classification errors, Fig. [22] reproduces 
the graph using only MK spectral types from Keenan & 
McNeil (1989) or Buscombe (1998), when available. A con- 
siderable scatter remains. Some of it is due to the real 
scatter in the properties of the stars (surface abundances, 
gravity, unknown variability). For supergiants in particu- 
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lar, and especially at low temperatures, the scatter also re- 
flects the large intrinsic uncertainties associated with the 
relatively poor quality of the model fits. We expect the 
spread to shrink once models with a wider range of pa- 
rameters (surface abundances, micro-turbulent velocities) 
will have been computed. 

We have also examined the diagrams of estimated T eB 
vs. spectral type obtained when any available optical data 
is included in the fitting procedure. They are similar to 
those described above. Individual stars are moved by up 
to 200 K, but no systematic trend can be clearly identified 
(because the stars that move most are also those for which 
the fits are poorest). Despite the added difficulty of fitting 
a broader wavelength range, the final dispersion is not 
significantly enhanced. 

7. Conclusions 

We have presented two grids of PHOENIX models for the 
spectra of luminous cool stars, one at solar metallicity, 
the other with RSG-spccific surface abundances. We have 
described the properties of these models and compared 
them with observations, with a focus on the molecular 
features found in near-IR spectra at resolution A/ A A ~ 
1000. At these wavelengths, red giants and supergiants 
dominate the integrated light of stellar populations. Our 
main conclusions are the following. 

— Models must be computed with a wavelength sampling 
step of about 0.1 A in order to reproduce the low res- 
olution near-IR spectra adequately. 

— The solar metallicity models provide a very good rep- 
resentation of empirical spectra of giants of class III 
and of a large fraction of the luminous giants of class 
II. As expected, RSG-specific abundances are found 
inadequate for the modelling of the bulk of the giant 
stars (they are rejected because they provide poorer 
fits and lead to a zone of avoidance in the derived Toff- 
distribution). RSG-specific abundances are favoured 
for some class II giants, which may have suffered mix- 
ing in excess of standard first dredge-up. 

— Red supergiant spectra of spectral types G and K, 
and of luminosity class lb (sometimes also lab) can 
be reproduced reasonably well. Serious disagreements 
remain in the case of very luminous (la and some lab) 
and of very cool supergiants (type M). RSG-specific 
abundances tend to improve the fits to strong CN 
bands, although the global effect on the fit quality is 
not as spectacular as one might have hoped. However, 
changing the surface abundance ratios has a significant 
impact on the derived effective temperatures (the ef- 
fect is larger than that found when moving from 0.5 Z Q 
to Zq). Therefore, it will remain necessary to account 
for this effect of stellar evolution in future model grids. 

— While it is easy (relatively) to produce good fits to the 
spectra of either the J, or the H, or the K band spectra 
of luminous cool stars, it remains more difficult to re- 
produce all their optical and near-IR molecular bands 



simultaneously. As a result, estimated stellar param- 
eters (T eff , log(g), Ay) depend on the spectral range 
of the analysis. The effects of changes in the surface 
abundances on these parameters also depend on the 
wavelengths under study. 
— The T eff scales derived from the comparison of a collec- 
tion of near-IR stellar spectra (1 — 2.4 fim) with mod- 
els are generally consistent with previous scales, albeit 
with considerable scatter. For cool red supergiants, the 
current uncertainties on individual estimated T eff val- 
ues frequently exceed ±100 K. 

About 20 % of the analysed red supergiant spectra have 
such strong CN bands that they call for models with 
high micro-turbulent velocities, and/or even more sur- 
face nitrogen than we have considered, and/or for grav- 
ities lower than \og(g) = — 1. The coolest of these are vari- 
able, and variability may contribute to the building of an 
extended atmosphere with low effective gravities. Large 
micro-turbulent velocities have been derived for a number 
of red supergiants in the past, and our first calculations 
confirm that increasing this parameter will help reproduc- 
ing the spectra of type la supergiants. In particular, a bet- 
ter agreement with observations is expected for the ratio 
between the first and second overtone CO band strengths. 
A grid of models is currently being calculated. Somewhat 
higher nitrogen abundances than we have explored are ex- 
pected to exist in nature, for instance when stellar rotation 
increases internal mixing. Because low resolution near-IR 
spectra of red supergiants are relatively easy to acquire, 
their comparison with models at the specific abundances 
predicted by stellar tracks with rotation will provide inter- 
esting tests of stellar evolution theory. Considering stars 
with lower gravities is a more challenging modelling task, 
as they will develop strong winds. In addition, the winds 
may be loaded with dust. Since winds are a well known 
empirical property of many red supergiants, the develop- 
ment of models that include winds is a necessity of the 
future. 
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